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Abstract 

Neutrino-nucleus cross sections for the detection of atmospheric neutrinos 
are calculated in relativistic impulse and random phase approximations. Pion 
production is estimated via inclusive A- hole nuclear excitations. The number 
of pion events can confuse the identification of the neutrino flavor. As a further 
source of pions we calculate coherent pion production (where the nucleus 
remains in the ground state) . We examine how these nuclear structure effects 
influence atmospheric neutrino experiments and study possible improvements 
in the present detector simulations. 
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I. INTRODUCTION 



Aside from the puzzling solar neutrino data there exist extensive measurements of atmo- 
spheric neutrinos, which are produced in cosmic ray interactions with the earth's atmosphere. 
Cosmic rays, mostly protons or a particles, hitting the atmosphere produce pions, which 
undergo the following decays: 



From this, one expects the neutrino flavor ratio. 



to be about 2. However, two existing experiments, 1MB and Kamiokande [Q], consistently 
measured a ratio of about 1 by studying charged-current neutrino-nucleus reactions [(i^e, e) 
and (z/^, /i)] in large underground water detectors. This could mean either there is a depletion 
of muon-type neutrinos or an enhancement of electron-type neutrinos. While the separate 
and z/e neutrino fluxes differ in various calculations, their ratio, however, seems to be largely 
model-independent , and the theoretical results are in clear contradiction with data. This 
discrepancy is known as the "atmospheric neutrino anomaly" . 

One explanation of this result could be the existence of neutrino oscillations of the muon 
neutrino into some other neutrino species 0] which would reduce the flux. However, 
before drawing this conclusion, one should also investigate more conventional sources for 
uncertainties in the measured flux ratio, which originate from nuclear physics and detector 
specifics. One crucial point is the experimental separation of electron and muon events. Most 
data are measured using water detectors. In a charged-current weak interaction with an ^^O 
nucleus, the incoming neutrino generates an outgoing lepton which is detected through its 
Cerenkov radiation. Electrons and muons are distinguished by the characteristics of their 
tracks. Electrons produce a showering "fuzzy" track, whereas muons are identified by their 
nonshowering track. 

Pions constitute a major background in this way of determining the neutrino flavors. 
In the experimental analysis, two- and more ring events, which signal the production of 
several charged particles, are rejected. However, a charged pion can be confused with a 
muon event if the lepton energy is below Cerenkov radiation threshold, and a neutral pion, 
with one decay photon missing, can be counted as an electron event. Therefore, a reasonable 
determination of the flavor ratio can only be achieved when the pion events are properly 
taken into account. 

One major source of pions is the decay of a A produced in the neutrino scattering. The 
delta can subsequently decay into a pion and a nucleon. An estimate of the number of pion 
events is only possible when the A-h response is properly calculated. Indeed, Monte Carlo 
simulation of the Kamiokande experiment estimated the number of pion events based 
on a rather complicated nonrelativistic model 0. However, atmospheric neutrinos have a 
broad energy spectrum ranging from a few hundred MeV to several GeV. It is therefore 
crucial to have a relativistic formalism that allows the calculation of the cross section for 
arbitrary kinematics. 
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As another source of pions, we consider coherent pion production and its role in the 
detection of atmospheric neutrinos. Coherent pions are produced in both charged- and 
neutral-current neutrino-nucleus interactions which leave the nucleus in its ground state. 
Traditionally, the coherent vr^'s are distinguished from the resonantly produced pions due 
to their strongly forward peaked angular distribution. However, in atmospheric neutrino 
experiments, due to a lack of directional information of the incoming neutrinos, it is not 
possible to distinguish the coherent tt^ events from incoherent pions. 

In the case of coherent charged pions the forward angle dominance of the cross section 
leads to a small opening angle between two charged particles {fi or e and charged pion) in the 
final state, and these two particles could be detected as one isolated electromagnetic shower. 
In contrast to incoherent charged pions whose Cerenkov tracks are usually identified as muon 
events, the special characteristics of coherent pions could lead to an increase of identified 
events and thus to underestimating the to ratio. 

This paper is organized as follows. In Sec. |I| we present the formalism for neutrino- 
nucleus cross section including p-h and A-h nuclear excitations within the framework of a 
relativistic mean-field model of the nucleus. Sec. ^ contains the formalism for coherent 
charged and neutral pion production. The resulting cross sections and production rates are 
shown in Sec. We close with a discussion and outlook in Sec. |V]. 



II. INCLUSIVE CROSS SECTION 



We outline the formalism for inclusive neutrino-nucleus scattering cross section including 
p-h (nucleon particle-hole) and A-h (delta-hole) nuclear excitations. We consider a neutrino 
with four- momentum k = {E^, k) which scatters from a nucleus via boson exchange 
producing a charged lepton with four- momentum k' = (£^k',k'). Using impulse approxima- 
tion, the double differential scattering cross section from a target nucleus with mass number 
A is given by (we assume a symmetric N = Z nucleus): 

_jPa_ _ _AG^co^W [^^^ n-] . (2) 



Here p = 2kp/37i'^ is the baryon density with Fermi momentum kp, 9c the Cabbibo angle 
(cos^6'c = 0.95), Gp is the Fermi constant. The leptonic tensor L^^j is defined as 

L^u = 8 {k^jkj^ + kyk^ — k ■ k g^^ =p i^ap^ivk k , (3) 

with the minus (plus) sign denoting neutrino (anti-neutrino) scattering. 11^ is the polariza- 
tion tensor of the target nucleus for the charged weak current. Here we consider p-h, Hp^, 
and A-h, H^"^, contributions to the polarization: 

nr = + nr. • (4) 

In the impulse approximation, the p-h polarization takes a simple form. 
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with the weak-interaction vertex given in terms of single-nucleon form factors parameterized 
from on-shell data, 

iq' = qo - q') (6) 

The form factors Fi, F2, Ga and Fp are given in the appendix of Ref. []^. The pseudoscalar 
form factor Fp is constructed from PCAC, and its contribution is suppressed by the small 
lepton mass. 

In an approximation where the nucleus is in a mean-field ground state of the Walecka 
model, the nucleon propagator G{p) is given by [§] 



G*{p) = {f + M* 



1 77T 



(7) 



where the effective mass M* and energy E* are shifted from their free-space value by the 
scalar (S) and timelike component (V) of the mean fields. 



M* = M + S] = + M*2 ; p*^' = (pO -V,p) . (8) 



In the Fermi gas approximation of non-interacting particles, V and S are set to zero. Analytic 
expressions for the imaginary part of the polarization [Eq. (^] can be found in Ref. 0. 

The impulse approximation can be improved by including RPA effects. An RPA calcu- 
lation uses the same expression for the cross section as in Eq. (|]) with the replacement: 

U^^^p^ = + AU^^^p^ . (9) 

An^p^ represents many-body correlations mediated by isovector particles, vr and p incor- 
porated with the Landau-Migdal parameter g' . 

The cross section involving A-h excitations is calculated similarly in Hartree approxima- 
tion using the A-h polarization n^"^. The weak interaction contains vector current (v) 
and axial-vector current (a) contributions, therefore, we split fl^"^ into 

nr. = i^Ahr + i^zr + i^AHT" + i^zr ■ m 

Here (n.^/^h)^'^ and (H^^^)^'^ are interference terms of the vector and axial-vector currents. 
These polarization tensors are written in terms of the spin 3/2 delta propagator and appro- 
priate weak vertices as, 

(nX,),. = - tj ^Tr[r^^(-g,-p) S^-ip) rUq,p) G{p-q) ] 

+ {hj) = {a,v). (11) 

The trace is taken for the Dirac matrices as well as the isospin matrices. For S'^'^{t) we take 
the Rarita-Schwinger form of the free spin 3/2 propagator with momentum t : 



The vector part of the nucleon-delta vertex has been studied in the case of the 7A^A tran- 
sition ITol , 

= V2F^T^V^, (13) 

with the isospin raising or lowering operator . The vertex for the axial A^A transition 
is given by PT}^TB[] 

ri = -"-^GATH9,u-^) (14) 

with r^A = fnNA/ fnNN ~ 2. The form factors and Ga are defined in the appendix of 
Ref. i. 

In mean-field approximation, the propagation of a A is modified by background scalar 
and vector fields similarly as the nucleon propagator. This can be understood from a chiral 
soliton model where the delta is a rotational excitation of the nucleon. That is, in the 
medium, the delta is expected to be influenced by the same magnitude of the vector and 
scalar potentials as the nucleon's because a and u are isoscalar. The calculation of the cross 
section proceeds in the same way as in the case of the free A. 

To include the delta decay width in our calculation, we average the nuclear response over 
the delta mass with a Breit-Wigner distribution [|l^]. The averaged cross section follows as 



(Pa r , 2 Po' 



(Ml- /.2)2 + Miri ^'^^ 

integrating from threshold to infinity. However, the decay width of a delta in the nuclear 
medium is not well determined both theoretically and experimentally. The ttN decay channel 
of the delta is partially blocked by Pauli blocking but there is an additional spreading width. 
Thus, determining the in-medium decay width is a nontrivial problem. Here, we assume the 
in-medium width to be the same as the free width Ta = 115 MeV. More theoretical and 
experimental work has to be done regarding this question, however. 



III. COHERENT PION PRODUCTION 

Coherent pions are produced as decay products of virtual p-h or A-h excitations of 
a finite nucleus. More specifically, the momentum transfer from the incoming neutrino 
virtually excites the nucleus through p-h or A-h and the nucleus decays back to its ground 
state by emitting a pion. Here we derive the formalism for the A-h excitation as it dominates 
the p-h excitations of the nucleus. The corresponding scattering diagram is shown in Fig. |l|. 
Note that the coherent pion is on-shell even though the momentum transfer in the initial 
scattering is space-like. The missing momentum is provided by the recoil of the nucleus. 

In a local density approximation, the finite- nucleus polarization Hfn(q', q'] ^) with incom- 
ing and outgoing momenta q and q', respectively, is approximated in terms of the nuclear 
matter polarization Hnm(q'5 q'] 
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n^N(?, d\ a;) = V F(|q - q'l) Jl^^{q, q'- u) , (17) 

where F(|q — q'|) is the elastic form factor of the target nucleus. We obtain F(|q — q'|) 
by Fourier transforming the nuclear ground state density obtained from a self-consistent 
relativistic mean field calculation for finite nuclei The volume for a target nucleus with 
mass number A is given by 



V 



2k'p 



As shown in Fig. |I], the incoming neutrino interacts weakly with the nucleus exciting 
an intermediate A-h state. The vertex of the weak interaction contains vector and axial- 
vector parts while the ttNA vertex, on the other hand, involves only the axial-vector part 
contracted with a pion four-momentum q'. Since the coherent-pion production cross section 
is dominated by pion momenta q' parallel to the momentum transfer q, there is negligible 
contribution from the part of the polarization that mixes the axial- vector and vector current. 
We ignore this small contribution throughout this paper. Consequently, the cross section 
for coherent pion production is identical for neutrinos and antineutrinos. 

The double differential cross section for coherent pion production is 

|k'||q'| Gl /•! 



,^ , , dcos6a' F'^(cos6o'] 



-^oo|nNMi 



|2 



+2Loi Re[nNM(nNM)*] + -^ii|nNMr 



(18) 



where is the angle between q and q'. Here we observed that the integral in Eq. ( pTSD is 
dominated by small angles, which implies that contributions from transverse components 
of the polarization, H^m and H^j^, are suppressed. Energy conservation requires the pion 
energy E^^ to be equal to the energy transfer to the nucleus go- The pion three-momentum 
q' is obtained from the mass-shell condition 



|q'| = \/g2-m2 . (19) 

The antisymmetric part of the leptonic tensor L^,^ does not contribute to Eq. ([I8|) , as it 
is contracted with the vector axial-vector mixing term of the polarization. The nuclear 
matter polarization, Hj^j^, can be evaluated for coherent charged or neutral pion production 
by taking appropriate vertices. The calculation of the real and imaginary parts of the 
polarizations H^m and H^im is summarized in Appendix. 

In a relativistic mean field approximation, the properties of the nucleon and A are 
modified by strong scalar and vector fields in the nuclear medium. These effects can be 
incorporated by modifying the nucleon and delta propagators as we described in the inclusive 
calculation. Since the axial-vector vertex does not contain the energy of the A or nucleon, 
the contribution from the constant vector mean-fields is eliminated in the calculation of the 
polarization by a simple change of integration variable. As the A is unstable, we include a 
width Fa by using a complex mass 

Ma ^ = Ma - i rA/2 , (20) 



in the denominator of the delta propagator This prescription yields similar results as 



the usual Breit-Wigner folding. Note, we have not included pion distortions in Eq. 
These may somewhat reduce the cross section. The effect of distortions and a possible A-h 
spreading potential remain to be investigated. 
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IV. RESULTS 



In this section, we present numerical results of our calculation. We start with the inclusive 
cross sections and their role in atmospheric neutrinos followed by a discussion of the coherent 
pions. In the case of the inclusive cross section, we discuss the case of muon neutrinos, but 
the general features of the results hold for electron neutrinos as well. In addition to the 
relativistic Fermi gas description of the nucleus, the effects of mesonic mean fields are also 
considered in the mean-field approximation (MFA). The target nucleus is assumed to be 
^^O with a Fermi momentum kp = 225 MeV. In MFA, we use the same scalar and vector 
couplings for nucleon and delta ("universal coupling"). For kp = 225 MeV, the effective 
masses are 

Ml = 931MeV ; M*^ = 638MeV , (21) 

and the vector self-energy is = 239 MeV. 

In atmospheric neutrino experiments such as Kamiokande and 1MB the cross section 
of interest is da/dE]^', which is obtained from Eq. (Q) by another numerical integration over 
the scattering angle of the outgoing lepton. Fig. ^ shows the cross section da/dE]^/ for 
E^ = 1 GeV. Note that the cross sections from the p-h excitation are peaked at a large 
lepton energy, which corresponds to a small nuclear excitation energy. This is the region 
where nucleon correlations are important. The nucleon correlations are modeled using RPA 
which includes a {p + tt + g') residual isovector interaction. The RPA calculation]^ in the 
mean-field ground state is also shown in Fig. |[ The cross section for stable deltas vanishes 
around i?k' = 0.7 MeV, which gives the minimum energy transfer for delta production. The 
curves including the delta width, however, spread over to higher lepton energies because of 
the averaging process. 

The mean fields reduce the p-h and the A-h responses by 30% at the peaks. However, 
RPA effects reduce the peak by 50% at high muon energies. The reduction of the A-h 
response by mean fields is interesting in the Monte Carlo simulation of atmospheric neutrino 
experiments. The strength of the A-h cross section provides the number of pions produced 
from the A decays. These pions may cause an uncertainty in determining the lepton flavor 
in the experiment. Indeed, the Monte Carlo simulation of the Kamiokande detector |^ 
takes into account the number of pion events based on a nonrelativistic model by Fogli 
and Nardulli 0. Their nonrelativistic results are qualitatively reproduced in our free delta 
calculations. Thus mean-field effects are not included in the existing simulations and could 
provide significant systematic errors. 

To get a rough estimate of the reduction in the number of events, we fold da/dE]^/ with a 
simple model of the atmospheric neutrino flux. We approximately fit the atmospheric muon 
neutrino flux with the formula 

0,^(E,) = 220 E;'-' ^ ^ , (22) 

m^ sr GeV sec 

and assume 



Note that our RPA calculations do not include A-h or mixtures of p-h and A-h. 
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0Pm = ^'^M = 201^. = 20P, 



(23) 



Using this flux, we calculate the yield Y in units of events per Kton year of detector exposure 
and per 100 MeV energy bin, 

Y = 1.1938 X 10^° I ^,^XE.,)^dE,^ , (24) 

and plot the result in Fig. |^ (a) with respect to the muon energy for the p-h and A-h 
excitations. Note, the A-h calculations include the decay width of the delta of 115 MeV. 
However, in the case of the inclusive process the quantitative results are insensitive to the 
decay width. At low muon energies, the mean fields in the impulse approximation reduce 
the total p-h events from 33 to 28 while including RPA yields 22 events (in a 100 MeV bin 
after one Kton Year of exposure assuming 100 percent efficiency). Thus, nuclear structure 
effects reduce the quasielastic charged-current events substantially. Since the electron events 
are also reduced by a similar amount by RPA correlations, this reduction does not directly 
affect the ratio r of Eq. (|1|). However, it can change the energy dependence of the cross 
section. 

Also, as one can see from Fig. ^ (b), the A-h contribution is about 40% of the RPA 
calculation at its minimum and, as the muon energy increases, the A-h response becomes 
more important. The A-h response is even larger than the RPA response at high muon 
energies (the ratio is 1.06 at E^/ = 2 GeV). Note that the A-h response is reduced by 10% 
to 25% due to mean-fields. Incoherent pion events should be decreased correspondingly. 
Therefore, the reduction due to the A mean-fields combined with the event reduction from 
RPA should be properly taken into account in the simulation. 

Not all deltas decay to vrN in the nuclear medium. In fact, in the medium the vrN decay 



is partially suppressed because of Pauli blocking |T5[. However, the A in the medium has 



additional channels of pionless decay. The coincidence experiments of {p, n) reactions at 



the A resonance |jT6[ showed that a large portion of deltas in the medium decay without 
emitting pions. Furthermore, in theoretical calculations, a value of about 70 MeV has been 
used as the 2p-2h partial width ||T8|JT7[] . Thus, experimental and theoretical results indicate 
a substantial amount of pionless delta decay in the medium. This 2p-2h decay mode of 
the delta has not been directly included in present detector simulations. At best, its effects 
have been partially incorporated using a very crude model of pion absorption without any 
explicit reference to A production. 

The pionless decay of the delta and the reduction due to the mean-fields, suggest that 
the simulations may substantially overcount the total number of pions. At this moment, 
we do not know how these effects change the result for the atmospheric neutrino anomaly. 
But the results indicate that the Monte Carlo simulations in current atmospheric neutrino 
experiments can be improved to treat pion events more accurately. 

Now we present results for coherent pion production in neutrino scattering. Coherent 
pions may be important in atmospheric neutrino experiments because they could increase 
the electron-like events and therefore reduce the ratio r of Eq. (p. A coherent neutral 
pion produced in neutral-current scattering (i/, u') immediately decays into two photons 
and may be counted as an electron-like event. Likewise, a coherent charged pion can also 
be confused with an electron-like event when the two charged particles in the final state. 
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muon (or electron) and charged pion, move with small opening angle and make an isolated 
electromagnetic shower. Indeed small-angle scattering dominates the coherent cross section. 

The single differential cross section da/dE^^ is obtained by numerically integrating 
Eq. (|T^) over the scattering solid angle. Again, we use the free width F^, = 115 MeV 
in our calculation. Pauli-blocking and additional decay channels for the A in the medium 
will change this value (and generate a non-trivial momentum-dependence). However, as 
theoretical in-medium values are still very uncertain, we will not explore modifications of 
the width at this point. More work in this direction has to be done in the future. 

In the atmospheric neutrino experiments, both ^x- and e-neutrinos participate in pro- 
ducing coherent pions. Furthermore, in contrast to the inclusive scattering process where 
the cross sections of antineutrinos are strongly suppressed, neutrino and antineutrino con- 
tributions to coherent pion production are of the same size. Combining the events from 
/i-type antineutrinos as well as e-type neutrinos and antineutrinos could enhance coherent 
pion production substantially, and it is of interest to quantitatively compare the coherent 
pion events directly with the inclusive charged-current events through p-h nuclear excitation 
for a realistic spectrum of neutrinos. 

Using the e- neutrino flux {(py^ = 0pj, we calculate the coherent pion events and compare 
them with electron production in scattering. The neutral pion events per kiloton-year are 

Y(7r°) = 6 X 1.194 x 10^°! <j){E,J^{E,JdE,^ . (25) 

The factor 6 includes the contribution from antineutrinos as well as a factor of 2 from /i- 
neutrinos. Note, for the purposes of illustration, let us assume all coherent pion events result 
in electron-like tracks. The actual identification of these events depends on the details of a 
detector. Fig. ^ (a) shows the coherent vr" events with and without mean fields and A decay 
width. The electron events from are obtained from 

Y(e) = 1.194 X 10^°| ^(E,jJf-^(E,JdEy^ . (26) 

The differential cross section da/dE]^i is calculated in relativistic impulse and in random 
phase approximation [0 . Note that the horizontal axis in Fig. ^ (a) shows either the coherent 
pion or the electron energy. 

Figure ^ (a) also shows that at roughly E = 300 MeV the coherent pion events, neglecting 
decay width and mean fields, are comparable in size to electron events calculated in RPA. 
This is certainly interesting because the electron-like events can be increased by these pions 
and therefore might help to explain the small ratio r of atmospheric neutrinos at ii^ = 300 
MeV. However, inclusion of mean fields reduces the pion events by a factor of 2. Furthermore, 
the A decay width exacerbates the situation: the number of events including mean fields 
and A decay width (thin dashed line) are only 10% of the electron events calculated in RPA. 

Next, we consider the coherent charged pion events. In this case, the Cerenkov light 
detector records the total energy of pion and muon (or electron) assuming they are moving 
with small opening angle and make an isolated single track. Overall energy conservation 
imphes 

E,^ = E^ + E^ or E,^ = E^ + E, (27) 
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depending on the incoming neutrino type. 

The coherent charged pion events are calculated from 



Y(7r±) = 2 X 1.194 x lO^Vl^^.e) / -^^dE, 



+ 4 X 1.194 X 10^V(i?.J / ^dE^ ■ (28) 

The first term results from e-neutrinos (the factor of 2 includes contribution from antineu- 
trinos) and the second term originates from muon-type neutrinos (with a total flux factor 
of 4). 

Figure ^ (b) shows the coherent charged pion events and the electron events calculated 
in RPA. As in the case of vr^s neglecting mean field and decay width of the A results in 
cross sections for the coherent pions comparable or even larger than the quasielastic ones. 
However, again mean fields and the decay width substantially reduce the pion events: At 
E = 1 GeV, the full calculation amounts to only 23% of the quasielastic electron events, 
but at E = 2 GeV, it increases to 30%. Since the pion events become important as E 
increases, our findings provide an important systematic error to the interpretation of the 



recent Kamiokande experiment in the multi-GeV energy range [|T9 . 

The effect of neutral and charged coherent pions can be substantial, the actual numbers 
are very sensitive to the delta decay width, though. Coherent vr" production could be large at 
E = 300 MeV while coherent charged pions may become important for high visible energies. 
To solidify our results, however, a more elaborate study of the A properties in the nucleus 
has to be performed. 



V. SUMMARY 

In this work, we have calculated neutrino-nucleus cross sections using a relativistic for- 
malism to investigate possible improvements for the Monte Carlo simulation of Kamiokande 
experiment. We have found that RPA corrections reduce the charged-current neutrino events 
by up to 37% for particle-hole excitations. Because of this reduction, the A-h response could 
be important and is found to give significant corrections to quasi-elastic nucleon knock-out 
processes. We have found that the relativistic mean-fields reduce the A-h responses by 
10% to 24%, which correspondingly reduce the pion events in atmospheric neutrino exper- 
iments. Combined with the additional channel of the non-pionic decay of the delta in the 
medium, our findings suggest possible uncertainties in the present Monte Carlo simulations 
of detectors and the interpretation of the measured data. 

We have also discussed coherent pion production in neutrino-nucleus scattering. We 
have shown that coherent pions can produce a sizeable background in the atmospheric 
neutrino measurements. As neutral pions and coUinearly outgoing charged leptons and 
pions may mimic electron-type events in the detector, this effect can produce uncertainties 
in determining the flavor ratio of the incoming neutrinos. Coherent pions may provide 
important systematic errors for the recent Kamiokande experiment in the multi-GeV energy 
region [|T3 . In order to get a clearer quantitative estimate of the influence of coherent pions, 



improved calculations should be done and these should be included in detector simulations 
of atmospheric neutrino experiments 
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APPENDIX: 



Here we evaluate the nuclear matter polarization, (nNM)/^, for coherent pion production. 
The polarization is defined as 



knNM)A 



Tr 



ri, S^-ip') Tiiq') G\p) 



^KMq"; 3 = {'^^y). (Ai) 

where p' = p + q. The superscript j indicates the polarization for the charged or neutral pion 



production. The axial-vector (F^,^) and pion (F^) vertices can be found from Ref. |11-13 
The noninteracting nucleon propagator in Hartree approximation reads 



ITT 



G\p) = (^+M)—S{po - E^)e{kF - IpI) . 
The polarization W [q) involves only axial- vector vertex and is given by 



(A2) 



d'^p 



Tr 



T,pS^"{p')r^, G\p) 



(A3) 



where V^y = Qav — lalv/^- The off-shell term {^alv/'^) insures that the axial vector vertex 
satisfies the relation "^^T = 0. The coefficient functions C^{q^) are given as 



1^ cos^, GA{q^) r^A for 



J =71^ 



for j = TT^ . 



(A4) 



Here 6c is the Cabbibo angle and (j^(g^) is the axial- vector form factor for the p-h excitation. 
The value of the axial vector A^A transition strength r^A is somewhat model dependent and 
we choose an intermediate value 2 in our calculation. 
After some algebra, Eq. (|A3|) becomes 



KM = -c'( 



' dE, 

M 



dx 



T, 



+ {la -Qa) ■ 



(A5) 



Here T^u is the result of the Dirac trace. 



T, 



^Vi^VvV + ^VvP^J) +P-P Qf^vP -IQp-p p^p^ 
-3Ml{pf,p^ + p^p'^ - hp ■ pgf,u) 
+2MMA{9g^, Ml - g^,p'^ - %p'^p\ 



(A6) 



and Bp is the Fermi energy Ep = \jkp + M^, and the x = cos6' is the angle between q 
and p. As long as A is stable, the angular integration can be done analytically while the 
remaining integration over Ep is performed numerically. A similar calculation in a mean 
field approximation can be done by replacing M and by the effective masses. The 
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contribution from vector mean field is eliminated by a change of variable in the energy 
integration. We include the A decay width by replacing the A mass in the denominator of 
Eq. (^) as 

Ma ^ = Ma - ir/2 , (A7) 
and the double integrations are done numerically. 
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FIGURES 



FIG. 1. Feynmam diagram for coherent pion production through A-h excitations. In the 
charged-current reaction, k' is the four momentTim of the muon or electron, and q is the mo- 
mentum of the boson. In the neutral-current reaction, k' is the momentum of the scattered 
neutrino and q the momentum of the exchanged . 

FIG. 2. Differential cross section da/dE\^i as a function of muon energy E\^i. The thick solid 
and dashed curves are for the free A-h calculations, without and with decay width, respectively. 
The thick dot-dashed curve is the free p-h calculation. The corresponding thin curves are for MFA. 
The short dashed curve is the RPA p-h calculation. 

FIG. 3. (a) shows events per kiloton-year versus muon energy. The bold solid line represents the 
free A-h calculation with decay width 115 McV, and the bold dashed line is for A-h in MFA. The 
thin lines are the corresponding p-h results. Events calculated in RPA are given by the dot-dashed 
curve. In figure (b), the solid curve is the ratio of mean field to free Fermi gas calculations for p-h 
excitations and the dashed curve is the similar ratio of A-h events. The dots gives the ratio of 
RPA to free p-h and the dot-dashed curve is the ratio of free A-h over an RPA calculation of p-h 
yields. 

FIG. 4. (a) shows electron events calculated in the Fermi gas model (short-dashed curve) 
and RPA (dot-dashed) involving only p-h excitations versus electron energy. Also shown are the 
coherent -k^ events calculated for A's with (bold long-dashed) and without decay width (bold solid). 
The corresponding thin curves are the cases including mean fields, (b) shows the similar curves 
for coherent charged pions. In this case, the events are plotted versus visible energy as discussed 
in the text. 
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